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The strong-isospin violation in ?/;(3770) — > J/V'+tt" via intermediate D meson loops is investigated 
in an effective Lagrangian approach. In this process, there is only one D-meson loop contributing to 
the absorptive part, and the uncertainties due to the introduction of form factors can be minimized, 
"^v^' ' With the help of QCD spectral sum rules (QSSR), we extract the J/il^DD* form factor as an 

implement from the first principle of QCD. The DD*-k'' form factor can be well determined from the 
, experimental data for D — !> tvIu. The exploration of the dispersion relation suggests the dominance 

■ of the dispersive part via the intermediate D meson loops even below the open charm threshold. 

This investigation could provide further insights into the puzzling question on the mechanisms for 
■!/;(3770) -5> non-DD transitions. 
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I. INTRODUCTION 



The non-DD decays of ip{3770) have attracted a lot attention during the past decades. As ip{3770) is the first 
state above the open charm DD threshold, its decay was believed to be saturated by the DD channel via the Okubo- 
' Zweig-Iizuka (OZI) connected diagram. Such an anticipation was supported by early experimental data which showed 
that exclusive decays of ^/'(3770) — > non-DD were negligibly small. Theoretical calculations of the perturbative QCD 
' (pQCD) leading order contributions also suggested rather small non-DD branching ratios for -0(3770) [H-Q. 
^ Interestingly, recent studies of the ■0(3770) non-DD decays in experiment and theory have exposed unexpected 

00 results which complicated the situation. In experiment, the DD cross section measurement by the CLEO collaboration 
suggests that the non-DD branching ratio is consistent with zero with an upper limit of about 6.8% [7-9]. Rather 
contradicting the CLEO results, the BES collaboration finds much larger non-DD branching ratios of ~ 15% in the 
" I direct measurement of non-DD inclusive cross section [l^]- Recently a next-to- leading-order (NLO) nonrelativistic 
, QCD (NRQCD) calculation of the cc annihilation width for -0(3770) suggests that the higher order contributions 
' can account for about 5% of the 0(3770) non-DD decay branching ratios at most [ll|. In Refs. [H, [l^, it was 
proposed that the open-charm threshold effects via intermediate meson loops (IML) could serve as an important 
. . nonperturbative mechanism to produce sizable non-DD branching ratios. Note that -0(3770) is close to the DD open 
E> ' threshold. A natural conjecture is that the DD threshold would play an important role in its production and decay. 
This mechanism turns out to be successful in the explanation of the decay of -0(3770) — vector + pseudoscalar as one 
of the non-DD decay channel of -0(3770) [H, [ll. 
. During the past few years, there have been observations of a large number of heavy quarkonium states [l^ at the 
B-factories (Belle and BaBar) and electron storage-rings (CLEO). Some of those states have masses close to open 
thresholds and cannot be easily accommodated in the framework of potential quark models. For instance, the well- 
established X(3872) is located in the vicinity of D*D threshold and its mass as a 1++ state is nearly 100 MeV lower 
than the first radial excitation of Xci in potential models. Such observations, on the one hand, have raised serious 
questions on the constituent degrees of freedom within heavy quarkonia, and on the other hand, raised questions 
on the role played by the open decay thresholds via the IML as an important nonperturbative mechanism in the 
understanding of the properties of those newly observed states. Such a mechanism symbolizes a general dynamical 
feature in the charmonium mass region, thus should be explored broadly in various processes. 

To gain further insights into the underlying dynamics and understand better the properties of the IML, we are 
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motivated to study the decays of -0(3770) — ?► J/ip + r] and J/ip + n^. First, we note that the decay of ip{3770) — )■ J/ip + rj 
is one of few measured non-DD decay channels in experiment with i?i?(-0(377O) — J/tp + 77) = (9 ± 4) x 10"'' fl^. 
One can estimate the branching ratio of ip{3770) J/ijj + tt" via tj-tt^ mixing based on the leading-order chiral 
perturbation theory. The mixing intensity can be expressed as jl6| 

° 4 TO^ - (to„ + 771rf)/2' 

Using Dashen's theorem one obtains 

Taking the rj — rj' mixing into account, the mixing intensity is slightly enhanced [l8l | 

e — eoA/3cos(/), (3) 

where v^cos^ = 1.34 would be unity if cj) is the ideal mixing angle. With the 77- tt" mixing intensity in a range 
of 0.01 ~ 0.02, the branching ratio of ■0(3770) J/ijj + tt'^ from Ty-Tr" mixing is at most the order of 10^^. This 
result actually sets up a limit for the rj-Tr^ mixing contributions in ^(3770) — >■ J/ip + tt'^. In contrast, the Particle 
Data Group (PDG2010) [H] gives an experimental upper limit Bi?(V'(3770) J/^/j + tt") < 2.8 x 10"''. A recent 
investigation of -0' "> J/'fp + tt" [l^ [1^ based on a nonrelativistic effective field theory (NREFT) suggests that the 
strong-isospin violation via the IML is relatively enhanced by 1/v in comparison with the tree- level contribution 
where the pion is emitted directly from the charmonium through soft gluon exchanges, where v ~ 0.5 is the velocity 
of the intermediate charmed meson. Since ■0(3770) is close to the DD threshold, we expect that such a strong-isospin 
mechanism would also play a role. As a consequence, the IML mechanism may lead to a sizable branching ratio of 
■0(3770) — J/ip + ir'^ which might be significantly larger than that given by the 7y-7r° mixing. In an early study [21m23i 1. 
the absorptive contribution from the intermediate DD in the decay of ^0(3770) — J/tp + •7r'^(7y) was calculated with 
an exponential form factor determined by the characteristic mass scale in the exchange channel. It was also argued 
that the real part contribution from DD and other heavier D meson loops should cancel each other in order to obey 
the OZI rule successfully in J/0 decay. However, it is found 24] that the IML effects may still be important in 
the decay of charmonium close to the open charm threshold. This is because the quark-hadron duality turns out to 
have been broken locally. As a consequence, the decay of a charmonium state can still experience the open threshold 
effects significantly if its mass is close to the open threshold. Such a scenario may imply that the real parts of the 
exclusive 1/^3770) decays could not be neglected and could explain the observed sizeable non-DD branching ratios of 
V'(3770) 0. 

In this work, we shall apply an effective Lagrangian approach (ELA) to investigate the IML effects in 0(3770) — >■ 
J/ip + 7r°, and demonstrate that the IML transitions have dominant contributions to this isospin-violating decay 
channel. As an important improvement of this approach, we shall implement form factors from QCD spectral sum 
rules (QSSR) for the off-shell J/tpDD* coupling vertex, while the D*D'k^ form factor can be extracted from the 
semileptonic decay of D — >■ tt^Iv. We mention in advance that this elaborate treatment will allow a reliable estimate 
of the absorptive amplitude of 0(3770) — >■ J/tfj + ir'^. Meanwhile, we also explore the real part of the IML contributions 
with the help of dispersion relation, within which the effective threshold is determined by the experimental data of 
V'(3770) ^ J/ij + rj. 

This paper is organized as follows: In Sec. [ll] the ELA for the IML transitions is formulated. In Sec. IIIII the form 
factors from QCDSR and D meson semileptonic decays are investigated. The numerical results for the isospin-violating 
decay V'(3770) — > J/ip + tt'^ are presented in Sec. IIVI and a brief summary is given in Sec. |Vl 

II. ISOSPIN-VIOLATING DECAY OF ip{3770) J/ip + tt" VIA IML TRANSITIONS 

A. Absorptive part 

As illustrated by Fig. [H only the DD{D*) loop (the meson in the parenthesis denotes the exchanged particle 
between J/^ and w'^) contributes to the imaginary part. The effective Lagrangians are as follows [25.. ,26.]. 

C^,"DD = -ig^"DDyj"''Dl df, A , (4) 

C^DD' = ^ 9^ A - Dj Df) , (5) 

Cd'D. = -igD'Dp{D'd^P,,Df ~D*;d^'P,,D^^) , (6) 
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where the couphng g^poD* = /{fip^/Mjj-Mi:,) GcV ^. The same convention has also been adopted in Ref. [24 
The decay ampHtude via the DD{D*) loop is 



E 

Polarization 



{2g^"DDP3. ■ ei){-2g^DD-e^uapq2^2 Ps^s)! 7^ — P ' ^5) 



V2 



Pi - Mf, Pi - Pi 



Ml. 



-f{pI 



(7) 



At each vertex of the loop diagram, the off-shell effect or the finite size effect should be taken into account by 
introducing the form factor F{pf), which can be regarded as the extended version of the local couplings in the original 
effective Lagrangian. The form factor is also necessary for cutting off the ultraviolet divergence in the loop integrals. 

^/*{92,e2) 



*(3770)(gi,ei 




FIG. 1: The _D-loop diagram contributing to the absorptive part. 



Applying the Cutkosky rule, the discontinuity of the decay amplitude is 

'i'*P5 w„ „ f „uu , P5P5 



Disc.Mfi -i{~2ni)'^ J j^--^{2g^,,DD){P3 ■ ei)eq^e-p.^^pty l-g^''' + 



Ml, 



xSipl - Ml) X Sipl - All) X ^V>z.^-bi)^^;z..o(pi) 

with tq^t'p^^i = ^ai3')^iQ_2^*2^p1 ^ short notation. Then, there are two form factors depending only on the virtuality 
of D* left in our calculation. In practice, the product of these form factors can be parameterized as the product of 
the local couplings and an empirical form factor: 

F^,^dd'{pI)Fd'dApI) = -2g^^DD'^^^j^F{pl) , (9) 



where a dipolc form factor is adopted. 



A2 - M2 ^ 2 



f{pI) - ( A2_pf ) ' (10) 



with A = Md* + akqcD and Aqcd = 0.22 GeV. 

Further deduction gives the discontinuity of the decay amplitude 



Disc.Mfi = 2ig^"DD J 



d^P3 F^DD'{pl)FD*D7,«{pi) 



{27r)32Ep,2Ep, pI - Ml, 

X27r5(p3 +pI - qi)eq2<L-p3pP3 ■ ei\constraints , (H) 



where Ep^ = V\Pi\'^ + Mf, and constraints = (p^ >0,pl> 0,pl = Ml,pl = Ml,p5 = 92 - P3,P4 = P5 + p). 
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Extracting the antisymmetric Lorentz structure Cq^e^upCiu, we get the tensor ampUtude 

x2wS{pI +P4 - 5?)P3P3|constraTOts (12) 

From Lorentz invariance, this tensor structure can be decomposed into terms built out of the external momenta and 
metric tensor: 

Disc.M% = g>"'MA+p''p''MB+q^q2Mc + {rq2+p''q^)MD, (13) 

where only Ma will contribute to the final result when we contract the tensor amplitude with the extracted antisym- 
metric Lorentz structure. 

Contracting the tensor amplitude with the metric tensor, we obtain 



d^P3 F^dd'(j>1)Fd'd^o{pI) 



)32£p32£p, pI - Ml, 

x2-k5{pI+p\ - ql)pl\constraints 

_ „. f |P3pd|p3| sinMMy) F^nD'{pl)FD,D^o{pl) 



y.2n6{pl+pl-q'l)pl\ 



constraints 

(14) 



In the rest frame of ?/'(3770) and setting the direction of 5*2 as the z-axis, the dynamic variables can be expressed as 

h = \q2\ 



r V V 



n = \P: 



\ \ 2M^n 



Ml„-AMl 



2 2 
V ^ ql = 



^_Ml„+Ml-Ml 



2M^n 

3 Ml„-Ml+Ml 
w = V = — 

q1 = V + w = M^ii 

Pi = ^ (15) 
Then, the Lorentz invariant amplitude is 

MI F^dd,{x)Fd,d„o{x) 
G{x) 

(16) 



g^^Disc.M^f'^ = ig^i^noD x -^-t — ; r x / da; ■ 



where the propagator of D* is G(x) = + - 2[v{v + w)/2- hrx/2] - M^,. 
Similarly, we can get other three Lorentz invariant amplitudes: 

p^p^DiscMl"^ ^ igv'DD X 

dx ■ 



87r(u + w) 

w{v + w)/2 + hrx/2fF^DD' {x)Fd* {x) 



G{x) 

(17) 
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q2nq2uDisc.Mfl 



/•I 



X 



87r(u + w) 

[v{v + w)/2 - hrx/2fF^DD, {x)Fd.d^o{x) 
G{x) 



(18) 



Plj.q2vDisC.MYi = ig^"DD X 



87r(?; + w) 

[viv + w)/2 - hrx/2][wiv + w)/2 + hrx/2]F.^DD.{x)FD,D^o{x) 

(19) 



Solving these four equations simultaneously, we can get the complicated expression of the invariant amplitude A4a 

and the absorptive part of the decay amplitude Disc.Ai fi. The charge conjugate contribution gives the same result. 

Because of the mass of ■(/'(STTO) being above the charmed meson pair, the coupling constant g^"DD and the isospin 
difference may be difficult to get from theory because of the rescattering mechanism. So we will extract this coupling 
directly from the experimental data: 



B. Dispersive part 

In principle, all the meson loops of which the thresholds are above the V'(3770) mass would contribute to the 
dispersive part (i.e. the real part) of the transition amplitude. Because of the introduction of form factors in the loop 
integrals, some model dependence seems inevitable in the evaluation of the real part. Given that the imaginary part 
of the amplitude can be reliably determined as in the previous subsection, we shall apply the dispersion relation to 
obtain the real part of the decay amplitude. Taking the assumption that the spectral density can be approximated 
by the extrapolation M.a{M1„) — )• Ma{si), we have the unsubtracted dispersion relation: 



2M^a{s 



the 

(2M^+)2 Si -M2„ 



i^dsi+P 



2M'l{s 



'Ms, 



(21) 



where M'a'^ corresponds to the charged or neutral D meson loop's contribution, and the factor 2 in front of M.a{s{) 
refers to the charge conjugate contribution. Then the total decay width is 



r = 



h 



4n 



Re[Mjr?+(^ 



2M'l°* 



127rM|„ 



Re[M) 



Tot 12 



Ml{v-w + 



V + w 



2i 



IE 



2M^ 
2i 



2\2 



{v + w) — {vw + h ) 



(22) 



where A^J"* — -""a "t ■"''a- 

Two points should be stressed: one is the upper limit of the dispersive integral, and the other is the virtuality 

dependence of the coupling g^"DD- Generally speaking, the upper limit of the dispersive integral should be infinity 
from the mathematical viewpoint. But in practice, we only take a finite effective threshold th because the spectral 
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density is only an approximation. It is presented that the form factor of D + J /ip {virtual) — )• is harder than that 
oi D + piyirtual) — ?> D because D can "see" the size of smaher J/i/; [23] ■ We expect that the heavier ^/;(3770) also 
gives a harder g-,p"DD{si) form factor at large si so that in a limited Si region the si-dependence can be neglected. 
There is also literature ;13| to take this si-dependence into account by adding a suppression factor exp(— /|p3p) into 
the integrand of Eq. ([2T|) . where / is the square of the interaction length ^28|]. We will discuss both points in detail in 
the following numerical analysis. 



III. FORM FACTORS FOR THE OFF-SHELL VERTEX COUPLINGS 
A. QSSR reanalysis of the form factor Ft/jdd* {pi) 

Since the mass of J/ip is below the lowest threshold of open charm DD, it is not possible to measure the form factor 
FtliDD-'ip'i) in experiment directly. There has been a systematic investigation of the charmonium to open charmed 
meson form factors in the framework of QSSR [1^ H^l- As a crucial criterion of QSSR, the pole contribution should 
take a dominant part in the dispersion integral. To our surprise, it seems not possible to satisfy this condition with 
the parameters given in the literature. This stimulate us to reinvestigate the F^y)y)*(p\) with the improved QSSR 
approach, and crosscheck the result with finite energy sum rules (FESR). 

We shall be concerned with the three-point correlation function: 

r^.fe,P3) = j cl4xd4ye^^'--e-^(^'^-'^)-^(0|r{j3(a;)j2t(y)jit(o)} |o) , (23) 

where Jl = C7i,c, = q'^^c and = iq^^c denote the interpolating currents for the incoming J/ip{q2, ^2), incoming 
-D*(P5j £5) and outgoing D, respectively. Taking the advantage of the unique Lorentz structure for the VVP coupling, 
we can decompose F^i, simply as: 



((J2,P3) = Hql,pl,pl)ef,^af3q^p'^ , (24) 



where ps ^ P5 + q2- The above expression has an arbitrary sign compared with that the preceding section. Using a 
double dispersion relation, one can express the invariant amplitude as: 

Aw^ J {s - qi){u - p^^) 

For the _D*-meson off-shell, the spectral density can be obtained from the Cutkosky rule presented in the previous 
section. 

On the phenomenological side, the three-point correlation function can be approximated by the lowest resonance 
plus the "QCD continuum" contributions, where the latter come from the discontinuity of the QCD diagrams from a 
threshold: 

V^(yi^) = Md {hU) + A , (26) 

and smears the contributions of all higher resonance contributions. In this way, the phenomenological part of the 
three-point function reads: 

KP^'^ = f^fi'^T^il'rfT ,.2^ + "QCD continuum" (27) 

[pl-Ml,)[ql-Ml){pl-Ml) 

where 6c = M'j^Mu^M^jfufD* f-ti;/fTic, the form factor with virtual D* is defined as 

m( MT^Ut/K ^^ {D*{P5)\4^0){D{p,)\D*{P5)J/i^iq2)) 

{D{P3)\J^J\J/'4'{q2)) = 2 — rf2 ' (28) 

Pb ^"d* 

with 

{D{p3)\D*ip5)J/^iq2)) ^ F^DD'ipl)ee,e,p,p, , (29) 
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and the decay constants are normalized as: 



{D*{p,)\Jl^\0) = MD>fD>e; 

{0\j'\DiP3)) 
{J/i^{q2)\Jl\0) 



rric 



(30) 



Matching the two sides of correlation function, and performing the Borel transformation (Laplace SR), the lowest 
perturbative diagram gives 



1 Pj-Ml, 
47r2 Sr 



ds du p{u, s, t) 



with 



p{s,t,u) 



3mc 

1 

2mr 



4m^ ^UJ^ 
{s-Ml)T^-{u-Ml)T2r 



3X2 



(31) 



-st + ml{s + 2t) ± y^s{s - 4m2)(t - 



(32) 



Aus, X2 = u + t— s + 2m^, and ti^2 are the inverse squares of the corresponding Borel 
>■ and T2 —T' 0, we obtain the FESR. Here, we neglect the numerically small gluon 



where t = p\, A = (m + s — i)^ 
masses. Taking the limits ti 
condensate contribution [30l |. 

To the leading order approximation where the three-point correlation function is evaluated, it is consistent to extract 
the decay constants fjj and from the corresponding two-point functions at the lowest order other than the value 
extracted from the experiment directly, e. g. .f ^t, = 0.405 ± 0.015 GeV. The QCD expressions of the pseudoscalar and 
vector two-point functions are well known [3l[. We show our analysis for fjj and in Fig. [5] and Fig. [31 respectively. 
Stabilities in both the two-point sum rule variables t^,_d and variation of the continuum threshold A are observable. 
We show in Fig. U the A behavior of the ratio t^/tq which is rather stable, especially for rUc = 1.26 GeV. The 
obtained optimal ratios (stable with A) are: 



7^ _ I 0.30 for rric = 1.26 GeV, 
7J^~|o.28 for mc = 1.47 GeV, 

while the ad-hoc phenomenological choice used in the literature is: 

^ = ^= 0.364. 

The relations between the two-point parameters t^^d and the corresponding three-point parameters ri_2 are [32 



~2' 



T2 



2 



(33) 



(34) 



(35) 



As follows, we adopt rric = 1.26 GeV and A^ = A^ = A as inputs since it is difficult to find a global maximum 
with rric = 1-47 GeV in our three-point SR and the variation of r ratio is more stable with rric = 1-26 GeV [s^. To 
obtain more concrete information about the form factor, we consider a large virtuality interval < — j)| < 5 GeV^, 
which is the same as in Ref. [2^. As an illustration, we show in Fig. [5] the form factor F^dd-, (p|) at p\ = —3 GeV^ 
with both A > 0.4 GeV and ti/t2 = 0.3. For simplicity, Jd* — 0.24 GeV is the same as in Ref. 29]. The ratios of 
the pole contribution versus the whole dispersion integral are also depicted in Fig. [5] and parameterized as 
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FIG. 2: (color online), (a) The two-point SR of fo versus td with Ad — 0.8 GeV. The red dashed line is for rric = 1.26 GeV 
and the green solid line for rric = 1.47 GeV. (b) jo versus A_d with the minimum of td adopted. 
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FIG. 3: (color online), (a) The two-point SR of versus with A^ = 0.8 GeV. The red dashed line is for nic — 1.26 GeV 
and the green solid line for rric = 1.47 GeV. (b) versus A^ with the minimum of adopted. 



With the increase of A, the pole contributions will become larger. It is obvious to see that the pole contributions 
are less than one half at the maximum ti = 0.05 GeV even with A as large as 1 GeV. The situation will be worse with 
larger D* virtuality. This phenomenon seems to be a common problem for the form factors of charmonium to open 
charmed mesons. To avoid this difficulty of the SR criterion, we deduce the form factors from Laplace SR varying 
with different A, and also show the predictions from FESR in terms of A. In principle, these two SRs should give 



the same solution, which means that the result at the intersection point is the reliable one from QSSR, see Fig. 6(a 
The form factor from the above method is shown in Fig. |6(b)| and we use three different parameterizations to extenc 
the form factor to broader regions of the D* virtuality: 
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FIG. 4: (color online). Behavior of the ratio of SR variables t^/td versus A for two values of rric = 1.26 GeV (red dashed line) 
and 1.47 GeV (solid green line). 



F.n^■(P^-3GeV^m=1.26GeV) 



i.On R(p,'=-3GeV',m^=1.26GeV) 




0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 



(GeV"') 

(a) 



X, (GeV"') 

(b) 



FIG. 5: (color online), (a) The n windows of the form factor and (b) the contribution from the corresponding pole at 
p| — —3 GeV^, where rric = 1.26 GeV, ti/t2 = 0.3, and A as a variable. 
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The form factor obtained in Ref. *29] with fixed = 0.405 GeV, fo = 0.17 GeV, = 1.3 GeV, = Ad 
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0.5 GeV and 0.09 < ri < 0.14 GeV ^ are parameterized by the Gaussian formula [2S 

-(-pi + 27)2" 



F^DD'ipi) = 19.9 exp 



345 



(40) 



Below the D* threshold, our improved form factors are slightly larger and decline slower than the one in Ref. |29| . 
In fact, the form factors used in our following calculation are usually restricted to a small region —5<pl < 2 GeV^ 
with on-shell D mesons. Therefore, those different parameterizations would not bring noticeable differences to the 
calculation results, although in a broader momentum region they turn out to be different from each other especially 
in the timelike region. Usually, pQCD predicts the power falloff of the form factors, we will use the dipole fit in the 
following calculations. Notice that there is no real roots, i.e. the unphysical state, in the denominator of our dipole 
fit (we label it as power fit hereafter to distinguish it from the empirical dipole form factor). 
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FIG. 6: (color online), (a) The intersection point of Laplace SR and FESR at pi = OGeV^. (b) The J/iIjDD* form factor 
derived from our method compared with the Gaussian fit !29j . The blue dash-dotted line is our monopole fit, the green solid 
line is our dipole fit, the red dotted one is our Gaussian fit, and the black dashed one is the Gaussian fit '2!||. The dots are 
our SR result. 



B. The Form factor of Fi:,*D^o{pl) 

As mentioned earlier, the form factor Fjj* jj^o (pi) is determined by the D meson semileptonic decays, i.e. D — > n'^liy. 
In the momentum transfer region < 2 GeV^, we expect that the D* pole has the dominant contribution. Thus, 
the transition matrix element can be expressed as 

{n°{p)\dj,c\D+{p,)) 
~ J2{^'>ip)D*+ip,)\D+ip,)){0\d^,c\D*+ip,)) , \ 

_ FD*D^o{pl)MD.+ fD.+ ( PbjjJPjlPb\ 

Pi- Ml,, \ MU ) ' ^ ^ 

where we have used the following definitions consistent with the effective Lagrangian: 

{n\p)D*+{p^)\D+{pi)) = Fn^DApl){p-4) , 

{0\djf,c\D*+{p5)) = MD'^fo'^e^^. . (42) 
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The transition matrix element of the weak decay can be defined as [3J| 

{7T'>{p)\d^^c\D+{p + p,)) = -^[{2p+p,)MpI)+pMpI)] , (43) 

where the form factor f+iPs) has been measured with high accuracy, and can be parameterized as a modified pole 
formula [3^1 

Compared with the part of the weak decay form factor definition, we obtain the needed form factor: 

2 V2,f+{0)MD. 
Fd'D7:o(P5) = 7 -r^ , (45) 



where ao = 0.21 for ao = 0.24 for D+ ^&\, and UiO) = 0.64 from the lattice QCD simulations [36|, l37| for our 
numerical calculation. Note that /+(0) from QSSR 13^ are consistent with the lattice result very well. The local 
coupling go'+D+Tr" can thus be extracted from the form factor at pg = M|), with //j. = 0.24GeV, i.e. 

go'Drroipl) = Fd-d^pI = A^D-) - 9.97, (46) 

which is slightly different from the value extracted from the decay of D* D + n, i.e. gD*D-K° = go* dp I ^/^ — 
17.9/V2 = 12.7 [11]. One can also extract this coupling from QSSR or QCD light-cone SR. However, both SRs suffer 
from their inherent uncertainties and the corresponding couplings from most SRs are nearly the same as what we 
derived from the weak decay form factor. One can refer to Ref. ^40] for a review on this issue. Another reason for the 
discrepancy of the coupling values is that in the momentum transfer region < pg < 3 GeV^ which corresponds to 
the experimental kinematics, the form factor may vary drastically near the pole position of D* . 

One could of course calculate the dispersive part with empirical form factors, but we must emphasize that in some 
diagrams containing more D* mesons, e.g. D*D{D*) loop in the vector charmonium decay to a VP final state, the 
empirical dipole form factor used by most of the references is not enough to suppress the ultraviolet divergence in 
the loops so that we need other more complicated form factors such as the Gaussian form factor. Thus, we leave 
the direct calculation with empirical form factor aside. We present the empirical dipole form factor (Eq. (jlOp) with 
different cutoff in Fig. [7] and compare it with our QCD-motivated form factors. It is obvious that our form factors 
favor a > 2, which is consistent with the value used in the study of X(3872) decays [4l| . 



IV. NUMERICAL RESULTS AND DISCUSSION 



The determination of the effective threshold th of the dispersive part is not a trivial task. If we know the full 
information of the spectral density, th should be extended to be infinity. Usually, th = {Mo + Md-Y is taken 
in the literatures as a natural cutoff (4l| on the assumption that the spectral density can be approximated by the 
extrapolation of the imaginary part. As an improvement, it is assumed that the corresponding effective threshold 
of ■0(3770) — > J/i) + rj containing u, d components should be the same as that of -0(3770) — ?► J/0 + 7r°. So, we can 
determine the effective threshold from the decay width of 0(3770) — >■ J/0 + r]. The flavor mixing scheme is taken 
in our calculation, and the mixing angle between \nn) = {\uu) + \dd))/\/2 and \ss) is ap = 9p + arctan\/2 ~ 38°, 
where 9p is the mixing angle between the flavor singlet and octet. Note that, as mentioned in the Introduction that 
the IML contributions are relatively enhanced by 1 /v in comparison with the tree- level contribution, it leads to the 
dominance of the IML contributions in -0(3770) — > J ftp + rj which is consistent with the study of Ref. jT2|. It can be 
understood that the isospin violation from the rj-ir'^ mixing is different from the IML transitions. In the latter, there 
is no rj pole contributions to the strong isospin violations. 

In the following calculation of the dispersive part, we will take the SU(3) flavor symmetry for the production of ss 
and qq component within the light pseudoscalar mesons. It means the dispersive part of ss is approximated by the 
average of the dispersive integrals of qq. 

The decay constants of g^"DD in our numerical simulation are listed in TableHJ where g^^lf^D^'^^ correspond to the 
lower bound, central value, and upper bound allowed by the experimental data p^ . From the experimental branching 
ratio Br(V'(3770) J/i) 77) = 9 ± 4 x lO-^flil, we obtain the corresponding 5th = th- {Md + Mo-f in Tables IH] 
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p/ (GeV^) 

FIG. 7: (color online). The QCD induced form factor F^DD*{p'i)Fj^-,j^^a{p%) (charged D) compared with the empirical dipole 
form factors. The red dotted line is our power fit, the black dashed line is the form factor obtained in Ref. [2!i], the green solid 
line, the orange dash-dot-dotted line and the blue dash-dotted line correspond to the empirical dipole form factor with a = 4, 
a = 2 and a = 1, respectively. 

and mil One can see that the frequently used natural cutoff 5th = is not supported by our calculation. Notice that 
we distinguish the charged and neutral D{D*) masses in our numerical calculation. To estimate the uncertainty from 
J/ipDD* form factor, we choose our power parametrization and the Gaussian form in Ref. [2^ for comparison. The 
si-dependence oi g^noD is also taken into account by adding the suppression factor exp(— J|p3|^) into the dispersive 
integral, where / — 0.4 GeV~^ is extracted from the charmonium mass shift [2^. It should be stressed that the 
dispersive integral with the suppression factor is not considered priority than the original one with a lower effective 
threshold from the phenomenological viewpoint. Moreover, both imaginary parts numerically decrease faster than 
1/si, so the unsubtracted dispersion relation used here is self-contained. 

In most of the parameter space, the dispersive part of the branching ratio is dominant over the absorptive part for 
J/ipri, while for J/^tt" both absorptive and dispersive parts are important. The difference between different J/tpDD* 
form factor is small, despite that the absorptive part of our power fit is systematically larger than the Gaussian fit in 
Ref. [29] as expected. As shown in Table |lll the branching ratios of J/ ipir'^ obtained from the corresponding effective 
thresholds are in good agreement with the experimental upper limit 2.8 x 10^^ [15]. As an interesting investigation, 
we take the threshold asymptotic to infinity with the suppression factor, and the corresponding branching ratios of 
J/ijjTT^ are still below the upper limit except for g^^Jfjjjj- In contrast, the asymptotic limits of J/i/jt] are far beyond 
the upper limit of the experiment 13 x 10^"'. Then, even taking the suppression factor into account, the spectral 
information from other resonances and continuum are still ambiguous, so that the asymptotic limit is questionable 
and the effective threshold is still necessary. 

It is essential to recognize that the isospin symmetry breaking with the vertex couplings is also an important 
dynamic source apart from effects caused by the mass differences between the charged and neutral D* mesons. In our 
formulation, the coupling g^"DD and form factor Fjjtjj^o{p^) are extracted from experimental data which suggest 
different values for the charged and neutral couplings, respectively. Since the form factor Fjj* jj^^o (p'^) can be better 
fixed by the experimental data, the results listed in Tables U |ll] and IIIII also reflect the effects from the isospin 
breakings of g^"DD- It is interesting to note that the larger absorptive contributions actually favor smaller difference 
between the charged and neutral g^"DD couplings, which is also observed in Refs. '55, f23i] considering the theoretical 
Coulomb correction for g^"DD- Taking into account the dispersive part, the central values of g^"DD give relatively 
small branching ratios for ip" — >■ J /ip -\- tt'', while deviations from the central values can produce larger branching 
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TABLE I: Different g^"DD from experimental data [HI ]. 





Br +D°) 


Br {tp" -^D+ + D-) 


g,!i"DODO 




Coupling Fraction (N/C) 


low 


0.47 


0.45 


12.26 


14.62 


0.84 


cen 
yil>"DD 


0.52 


0.41 


12.90 


13.95 


0.92 


oT, 

y^" DD 


0.57 


0.37 


13.50 


13.25 


1.02 



TABLE IL The branching ratio without suppression factor. "Abs" denotes the absorptive part, and "Tot" is for the sum of 
the absorptive and dispersive part. The flavor mixing angle is ap = 38°. 





Our power fit 


Gaussian fit [29] 


low 

y^"DD 


cen 

yi>"DD 


iljlj" DD 


low 

yi>"DD 


cen 


y-ilj" DD 


Br^(Exp)(xlO"*) 


(5, 13) 


(5, 13) 


(5, 13) 


(5, 13) 


(5, 13) 


(5, 13) 


Br^(Abs)(xlO-'') 


2.13 


2.16 


2.18 


1.12 


1.14 


1.15 


(5t/i(GeV2) 


(0.17, 1.25) 


(0.17, 1.25) 


(0.17, 1.25) 


(0.9,2.5) 


(0.9,2.5) 


(0.9,2.5) 


Br^(Abs) (xlO-^) 


0.145 


0.320 


0.558 


0.074 


0.164 


0.287 


Br^(Tot)(xlO~*) 


(0.46, 1.36) 


(0.33,0.35) 


(0.71, 1.13) 


(0.49, 1.60) 


(0.17,0.23) 


(0.50,0.81) 



ratios for -0" — > J jil) ^ ifi . Within the present experimental uncertainty bounds [15], the predicted branching ratios 
for i/)" J /ip + 7r° are at the order of 10~^ ^ 10^"*. Confirmation of this decay branching ratio would be a strong 
evidence for the open charm threshold effects in ip" — > J/tp + ir'^. Note that our prediction of the absorptive part is also 
close to the prediction of Ref. ^ with isospin / = for -0(3770), i.e. Br^(Abs) = 8 x 10"^ and Sr^(Abs) = 2 x 10"^ 
which is a consequence of the similar values of the form factors in both approaches. 

V. CONCLUSION 

The ELA is very useful to investigate the nature of the near threshold charmonia and charmoniumlike resonances. 
The largest uncertainty of the ELA comes from the determination of the off-shell effect, i.e. the form factors. In 
this paper, we investigate the isospin violating decay of 0(3770) — J/0 + 7r°. In this process, there is only one 
D-meson loop contributing to the absorptive part, and the form factors encountered in the loop calculation can be 
relatively well controlled. With the help of QSSR, we extract the J/tfjDD* form factor as an implement from the 
first principle of QCD. The DD*tt° form factor can be well determined from the experimental data of D ^ ttIv, 
which has been measured with high accuracy. We also explore the dispersion relation to evaluate the dispersive part 
of "0(3770) non-DD decays, and find they take an important part in most of the parameter space. It means that the 
IML effects below the open charmed meson threshold cannot be neglected in general. Different from the traditional 
natural cutoff of the effective threshold in the dispersive integral, we extract them from the experimental data of 
■0(3770) — >■ J/ip + T]. Our prediction of the branching ratio of 0(3770) — )■ J/ip + tt" is less than 3 x 10~^ with the 
couplings g^i> z) d extracted from the central values of the data. Within the experimental uncertainty bounds for the 
extracted g^/zjjo^o and g,p"D+D-: the branching ratio of 01(3770) — )■ J/0 + 7r° can reach the order of 10~^. Notice 
that the understanding of the isospin violation of 5^// 00 is not a trivial task despite the Coulomb correction favors 



TABLE IIL The branching ratio with suppression factor. "Tot" is for the sum of the absorptive and dispersive part, and 
"Asym" means we take the asymptotic limit 5th 00. Flavor mixing angle is 38°. 





Our power fit 


Gaussian fit [29] 


low 


cen 
9ip"DD 


yip"DD 


low 

yi,"DD 


cen 

yiij"DD 


up 
ilij" DD 


Br^(Exp)(xlO-'') 


(5, 13) 


(5, 13) 


(5, 13) 


(5,13) 


(5, 13) 


(5, 13) 


Sth{GeY'^) 


(0.26,1.60) 


(0.26,1.60) 


(0.26,1.60) 


(1.15,3.50) 


(1.15,3.50) 


(1.15,3.50) 


Br^(Asym) (xlO^*) 


84.7 


84.7 


84.2 


35.2 


35.2 


35.1 


Br-^(Abs) (xlO''') 


0.145 


0.320 


0.558 


0.074 


0.164 


0.287 


Br^(Tot) (xlO"'') 


(0.48, 1.49) 


(0.33,0.37) 


(0.69, 1.05) 


(0.53, 1.73) 


(0.18,0.26) 


(0.47, 0.73) 


Br^(Asym) (xlO"*) 


13.76 


1.47 


2.86 


5.58 


0.62 


1.25 
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the experimental central value [42|. It is also suggested that a small admixture of isovector four-quark component 



of V'(3770) may also give a measurable decay rate to J/tpir^ Q- In Ref. [J| the four-quark component is viewed as 
a reannihilation effect of DD. To some extent, the nature of ■(/'(STTO) hides in the coupling g^"DD- Meanwhile, the 
forthcoming BESIII measurement of "0(3770) — J/?A -I- 7r° will be able to provide useful information about the QCD 
motivated form factors and clarify the role played by the IML. We plan to discuss the isospin violations with the 
charged and neutral couplings g^"D+D- aiid g^nj^oj^o elsewhere. 
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